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A hydromagneto-ionic theory has been developed within the framework of the Burgers 
formalism | L958a] which is a microscopic theory based on moments of the 4 Boltzmann trans- 
port equation. The effects on electromagnetic wave propagation of electron-electron and 
electron-proton collisions have been considered to the order of Chapman and Cowling's 
second approximation. The present theory is an extension of the magneto-ionic theory 
derived earlier by one of the authors [Wyller, 1961] for a fully ionized hydrogen plasma, in 
that it includes the effects of ion motions. The theory is therefore applicable to the high as 
well as the low frequency modes of wave propagation, such as the whistler, Alfven, and 
the retarded magneto-acoustic modes. Expressions have been derived for the refractive 
index, absorptivity, the wave polarization, and the zeros and infinities of the refractive 
index. Numerical applications have been given for the four characteristic modes of low 
frequency wave propagation, viz, the whistler mode, lower hybrid frequency, ion gyro- 
resonance, and the hydromagnetic mode. Applications of the theory to the solar corona, 
and future extensions to the terrestrial ionosphere have been indicated. 



1. Introduction 

The present paper considers the effect of the ion 
motions on electromagnetic wave propagation in a 
uniform plasma with a superposed homogeneous 
magnetic field. The conventional magneto-ionic 
theory considers the motions of the electrons only, 
which is a valid approximation for the high wave 
frequencies used in ionospheric propagation or for a 
'weakly ionized gas. When the wave frequency goes 
below the electron gyro frequency, the extraordinary 
electromagnetic wave becomes the oblique Alfven 
(whistler) mode, and for a further decrease below the 
ion gyrofrequency the ordinary electromagnetic 
wave becomes the retarded magneto-acoustic mode. 
Denisse and Delcroix have given an excellent illustra- 
tion of this in their book [1961a]. We note that if 
we take account of the pressure tensor, additional 
modes will appear. We have not considered these 
j nodes in our paper. 

Our work differs from most o treatments [Schluter, 
1950, 1951; Dungey, 1951; Astrom, 1951; Hines, 
1953; Piddington, 1956; Gershman, Ginzburg, and 



1 Presently at the Institute of Theoretical Astrophysics, University of Oslo, 
Norway. 



Denisov, 1957; Lehnert, 1959; Fejer, 1960; Denisse 
and Delcroix, 1961b; Allis and Delcroix, 1963] in 
that it takes account of collisional damping by a 
microscopic kinetic theory with a proper averaging 
over the velocity distribution function, of like 
particle interactions (proton-proton and electron- 
electron collisions), and of the heat-flow equations. 

We have seen that a microscopic treatment through 
the Boltzmann equation led to significant departures 
in the propagation factors for electromagnetic waves 
in the terrestrial ionosphere [Sen and Wyller, 1960; 
Leinbach, 1962]. Later on, one of us [Wyller, 1961] 
applied the Burgers formalism [1958a] to find the 
electrical conductivity for a completely ionized 
hydrogen plasma. The Burgers formalism proceeds 
via moments of the Boltzmann transport equation 
and the associated collision integral. It includes the 
effects of like particle interactions and of the electron 
and ion heat flows, and is equivalent to the Chapman 
and Cowling second approximation [1958]. The 
present work is an extension of Wyller' s [1961] work 
to include the effects of ion motions. 

In a recent work Kantor [1963] has considered the 
effects of ion motions on magneto-ionic theory with- 
out collisions and considered the properties of 
different directions of wave propagation with respect 
to the magnetic field. Our work checks with 
Kantor's formulae for the case of no collisions. 
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2. Basic Equations 2 

It is well known that the presence of a magnetic 
field in a plasma introduces an anisotropy in the 
propagation of electromagnetic waves through the 
medium. This anisotropy is characterized by the 
tensorial representation of the associated physical 
quantities governing the modes of propagation; 
namely, the electric conductivity and dielectric 
constant of the medium. After introducing the 
dielectric tensor (which is easily derived from the 
conductivity tensor) into the Maxwell equations, 
one can derive expressions for the refraction, 
absorption, and polarization of the propagating 
waves. With this end in view we shall set up in 
this section the fundamental equations that will 
enable us to compute the a-c conductivity tensor of 
the magneto-ionic medium which may transmit 
electromagnetic waves associated with the oscilla- 
tions of the ions as well as electrons. In doing so 
we shall adopt, as we have already said, Burgers' 
formalism [1958a]. 

Choosing a rectangular coordinate system (o, x, 

y, z) such that the 2-axis lies along the uniform 

-> 
applied magnetic field H , the electric conductivity 
tensor can be put into a very simple form. The 

electric field due to the propagating wave is rep- 

-> 
resented by an a-c field in that medium, E cos a>t = 
->-*-> ~? ~? -* 

(E x i-{-Eyj-{-EJc) cos cot, where i, j, and k are unit 

vectors along the x, y, and z axes and co is the excita- 
tion frequency of the wave, which is assumed to 
remain constant during the propagation of the wave 
inside the medium. 

We shall assume the plasma medium to be homo- 
geneous with no pressure and density gradients. 
We shall consider a fully ionized hydrogen plasma 
containing electrons and protons in equal number 
densities, ni=n 2 , throughout the medium. The 
subscripts 1 and 2 will be used consistently to 
denote protons and electrons respectively. 

The assumptions made above bring about a con- 
siderable simplification in the Burgers formalism. 
There is a simple relation between Burgers' "friction 
coefficient" Ki 2 =K 2 i=K and the mean electron- 
proton collision frequency v, namely, 



f ve- n *l /2kT vidv 2 
K=n 2 m 2 v, where v=- 



r 



e -m>\WTtfA v 



or 



2 co|e| ( m 2 y* _3JcT 



(1) 



(2) 



Here r D is the Debye length, and oo is the electron 
plasma frequency. 

In order to assess the full implications of the 
momentum and heat transfer equations of the elec- 

2 Sections 2, 3 and 4 were performed by one of us (A. A. W.) while on summer, 
leave in 1963 from the Institute of Theoretical Astrophysics, University of Oslo- 
Norway, to Williamson Development Company, Inc., West Concord, Mass. 



trons and protons, we state here some of the physical 
parameters of Burgers' theory that are relevant to 
the conditions of our problem. For a binary gas 
with known relative velocities of the colliding 
particles "g" the collision cross section is given by: 



S$=S$>=2* £*(l-oo* l x) 



bdb. 



(3) 



where b and X are Chapman and Cowling's [1958] 
collision parameters. In particular, 1=1 for momen- 
tum transfer with which we are primarily concerned 
here. The average collision cross section for all 
possible relative velocities is 



7(W - 



=Z2» = 



■yjir 



G*/2*T) J 



f 

Jo 



g2J+3 e -^WT S (l) dg . (4) 



Here the subscripts ' V and "t" refer to the species of 
the colliding particles and fi is the reduced mass. 
The superscript "Z" refers to the angular dependence 
while "j" refers to the velocity dependence of the 
collision cross section. In terms of these average 
collision cross sections, the quantity defined by 



-1- 



2Zf 



5 Z<}» 



(5) 



reduces to the value 3/5 for Coulomb interactions. 
In order to be able to estimate the effects of 
collisions between like and unlike particles sep- 
arately Burgers [1958b] introduces the following 
collision factors which turn out to be proportional 
to the average collision frequencies, 



m x m 2 



~5(m 1 + m 2 ) 2 Z 1 ^ 1) 



(electron-proton collisions) 



ki=- ( — ) ^ny (proton-proton collision.) (6) 

5 \ Zwfi 2 J Zji2 



K 2 ~- 



\ 2m 1 J 



12 y$ 



7(11) 



(electron-electron collisions). 



It will be seen later that the auxiliary quantity f 
which incorporates the effects of electron-proton 
and electron-electron collisions enters into the elec- 
trical conductivity tensor. It is defined by 



^ = k+k 2 +jY 2 ; y-- 



m 2 



(7) 



and Y 2 is a function of k only. 

Though the proton-proton collisions do enter into 
the equations of the heat flow through the factor 
«ci, their net effects in the electrical conductivity 
tensor vanish [Spitzer, 1956; Marshall, 1957]. On 
the other hand, in thermal conductivity the dissipa- 
tive effects due to ion-ion encounters is not negligible 
[Rosenbluth and Kaufman, 1958; Wyller, 1963a]. 

Burgers' momentum and heat flow equations in 
their time-dependent forms will be our fundamental 
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set of equations for the electrons and protons. We The gyrofrequencies are defined by 
will rewrite Burgers' equations [1958c] as follows. 



Momentum equations for electrons: 

ioyu2z— T 2 x+S2U 2y =v(u lx --U2x)+2pr2 X --zyvr 1X) 

ioiU 2y —Y 2y — s 2 u 2x =v{u ly —u 2v ) + zvr 2y — zyvr lyy 
ia)U 2z —T 2z =v(u lz —u 2z ) + zvr 2z —zyvr lz . 

Momentum equations for protons: 

ioM lx —T lx —Siii lv = —yv(u lx —u 2x ) — zyvr^-ir zy 2 vr ]x , 
iccUiy—T ly +s 1 u lx =—yv(u l y—u 2 y)--zyvr 2 y-\-zy 2 vr lyf 



rrii 



The foregoing set of 12 equations in 12 unknowns 
will now be solved to obtain the current flow compo- 
nents (um—u 2 h). Then the conductivity tensor 
will be found from the usual formula : 



(8) 



I=—n 2 e 2 w=<T'g, 



(14) 



3. Solution for Electrical Conductivity Along 
the Magnetic Field 

iuu u -T u =-yv(u u -v 2z )-zyvr 2z +zy 2 vr l2 . The momentum and heat flow equations for both 

(9) electrons and protons along the z axis are : 

Heat flow equations for electrons : . 

ioM 2z —T 2z =vtu l2 —u 2z ) + zvr 2z —zyvr l2j (15a) 

iuu l2 — T l2 =— yv{u u — u 2z ) — zyvr 2z +zy 2 vr lz , (15b) 



5 27 

io}r 2x +s 2 r 2 y=— ^ ^(u lx —u 2x )—ivr 2x +^ V>r lXy 

5 27 - 

ia)r 2y — s 2 r 2x = — ^ zv{u ly — u 2y )— frr 2y +— yvr ly , 



5 27 

iur 2z =— - zv(u lz —u 2z ) — &r 2z -\- — yvr lz . 



Heat flow equations for protons: 



(10) 



5 27 

iwr 2z = —- zv(u lz —i( 2z )—&r 2z +— yvr lz , (16a) 



5 _ 27 

iwr l2 ^zy 2 v(N lz —u 2z )—yK l ur 2z + Y {) y 2 ur 2z . (16b) 



According to Burgers we neglect the terms zyvr l2 , 

27 
zy 2 vr l2 , and — 7^^ for yr u <^<^r 2l . This can be 

verified by solving for r ]h and r 2h from the basic 
equations of the last section. 
From (16a) we have 



„ _ 5 ^ (u lz —u 2z ) 
2 (^+2w) 

Upon inserting (17) into (15a, b) we find that 



(17) 



-/ \ 5z 2 v 2 (ui 2 —u 2z ) , i0 N 

%uu 2 — T 2 z=v(u u —u 2z ) — o />- i ■ x — (18a) 



5 27 _ 

i«riif+ Vi«=2 ^{^\y—^2y)—y^r ly -\-— y 2 vr 2y , 

5 27 

i«r u =2 £7 2 K^iz— ^2z)— yki^u+Jq 7^22. 

(11) 

In the preceding equations the quantity ^co appears 
because of the time dependence of the particle 

velocity and heat flow velocity respectively induced 

-» -> 

by the alternating electric field §=& e iut . The 
mean flow velocity component, u sh , is the mean 
value of the particle velocities, % sh , 

u S h=£ S h' Defining the current flow velocity as w=u l —u 2 

(14) we have 
Also the accelerations due to the electric field are w 2 =^i 2 — 1£ 22 . (19) 

denoted by 

Thus we get from (18a, b), 



ioou l2 — T l2 = —yv(u lz —u 2z ) + 



2 (&> + i«) 

572 2 ? 2 (^i2-^22) 



2 (&+i<a) 



(18b) 



Ti== ^£> e i= + e '> e *=- e ( 12 ) i^W z +(l + y)T 2z =-(l + y)vW 2 + (l + y) £ 



5 z 2 Pw z 



so that 



ri=-yr 2 ,7= 



m 2 
mi 



(13) 



2 (&+io>) 

(20) 

Now solving for w z with separation of real and 
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imaginary parts we obtain 

+ tco(l + (l+ 7 ) |^~- 2 )]= = -(l+7)r 22 . (21) 

At this stage we introduce simplifying notations to 
condense the final expression for the conductivity 
tensor component along the z -direction and to eluci- 
date therein the effects of electron-electron and 
electron-proton collisions; we set 



8,=t: 



5 z 2 ?P 



2 £ £ 2 ? 2 +o> 2 



and 



V ^ 



(1— «b)5 



[l + (l+7)|] 



(22) 



The component of the current vector in the s-di- 
rection is then s;iven bv 



I z =—n 2 e 2 w z 



(l+7)q- n g, 



[l + (l+7)|][(l+7)i'3+i«] 



*o=£ (23) 

and the conductivity tensor component in the 2-di- 
rection is finally obtained 



U+tVo 



[i+(i+t) |] [(l-frK+ic 



(24) 



When we compare this tensor component with 
the corresponding Lorentz tensor component <r Q / 
(?+iw), we find that the ion-motions affect the con- 
ductivity along the ^-direction through the factor 
(1+t) while the electron-proton and electron-elec- 
tron collisions through 5 3 and vs. In particular, the 
expression (24) checks with Spitzer's d-c conduc- 
tivity [1962], when y and co are set equal to zero. 

4. Solutions for Conductivities Perpendicular 
to the Magnetic Field 

The magnetic field couples the x and y components 
of the momentum and heat flow equations, making 
the solution for the flow components u ix and u iv 
somewhat awkward to obtain. We may avoid this 
difficulty by introducing a coordinate system de- 
fined by the unit bi-vectors [Menzel, 1961]: 



<x=-i= (i+ij), 0=j= (i—ij), J=k. 



Now we may express the electric field as, omitting 
the time dependence for simplicity, 



-> -> -» -> | 

(b = fe x^~T fe yj\ Q> zk—~~= ( (b j 

V2 



(8x+ig v )P+gzy 9 



gia + &rP+g zJ. 



(26) 



Allis used a similar coordinate system in his for- 
mulation of the conductivity tensor for plasmas 
[Allis, 1956]. Si corresponds to the left rotating field 
component and & r to the right rotating component. 

We also write the flow components in this co- 
ordinate system 

•\2 y2 

We begin by solving for w r from the rewritten 
momentum and heat flow equations 

ic0U 2 r— r 2r — is 2 U 2 r = v(Ui r — U 2r ) + Zvr 2r (28a) 

icou u — Ti r +isiUi r = —yv(ui r —u 2r ) — zyvr 2r (28b) 



iur 2r —is 2 r 2r = —-zv(ui r —u 2r )—£vr 2r (29a) 

5 27 _ 

icor lr +is 1 r lr =-zy 2 p(u lr —u 2r )—yK 1 vr lr + --y 2 pr 2r . 



(29b) 



2 ~, 'V — 

From (28a) we have 

v(uir—u 2T ) = i<j)U2T—T 2r — is 2 u 2 r—zvr 2r , (30) 
which yields r 2r when put into (29a) 



r 2r 



i(s 2 —O^U 2 r + T 2 r 

'U-^z 2 )-i(s 2 -c 



(31) 



Putting expression (30) for v(ui r — u 2r ) into (28b) 
we obtain 



u lr =y 



(s 2 —oo) 
(«i+w) 



^2r 



(32) 



since r lr = — yT 2r . The electric field and heat flow 
(25) terms have dropped out. At this point we may 
return to (28a) with these expressions for u lT and 
r 2r and solve uniquely for 



u 2r — 



— [£?•— i(s 2 — ctf)]r 2 



{fe-co) 2 +[?(^-| S 2 )-i( S2 -co)][7^{|q^-?] + ^fe-co)} 



(33) 
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The total current flow associated with the right Now we introduce a new symbol 8 { and the expres- 

rotating electric field is then the difference of the sion for w r becomes 

electron and proton flows , 

{—a)T 2r 



Wr=- — j 7-T 

^=% r -^=[rg^-i]^- (34) -5(l-«+.-(— *)(l+^) 



For subsequent algebraic simplification we intro- 
duce the symbols 

a[£v— i(s 2 — w)]r 2 r 



«.=o 



?? 



2 H*?+ («-**)* 



(38) 



We can repeat the foregoing analysis and obtain 
the flow for the case of the left rotating electric field 



which gives 



Wj = 



-m 



[(s 2 — to) 2 — avA-\-i£v(s 2 — u>)] 
Expression (35) is equivalent to 

(+a)T 2r 



(35) 



|«(l-*a)+*(ft)+*2)(l+/»j) 



(39) 



where, as before, we have simplified the expression 
with the symbols 



w r ~- 



5z 2 



?v 2 



[(s 2 — co) 2 — avA+i£v(s 2 — co)] [&+i(s 2 —a))} 
[&+i(8 2 —w)] [&— i(s 2 — co)] 



, and 0=1+7 



(co+s 2 ) 



(36) 

On separating the real and imaginary parts in the 
denominator, we have 



° 2 2H 2 ? 2 +(co+s 2 ) 2 " 1 " i " * ' '(«-«i) 
From the foregoing, we obtain the total current 

/= — n 2 e 2 (wi a+w r @+ w z y). (40) 



w r = 



(+«)r* 



(37) 



Then the complete electrical conductivity tensor 



Rather than directly solve for u, and u we fa the hydromagheto-ionic theory becomes 
derive the electrical conductivity tensor in a lorm 
analogous to that of the Lorentz tensor form [Sen 
and Wyller, I960]. We first introduce the symbols 



7*= 



Wi 



w T 



Iz 



W z 



ri r -p i z -p 

21 — 1 2r — 1 2z 

which let us rewrite expression (40) as 



(41) 



_> (Tj+7r) i(ir,— 70 

*=SH -»(^-7,) (7.+7r) 








(45) 



where 



"Br 



^i{i + ij)+^yrr 2 S-ij)+y z T 2 Jcl (42) 
and 



(,+«§) 



[a^i + i(w — $2)] 



2 

w , 



7=/ a t+/ y i+i 8 fc=^[7i(g*-igy)i+T,(g,+igx)i 



7i= 



+7 r (gx+ig y )i+7r(g,-igx)i+27zg«fc], 
or, in matrix form 



(43) 






I, 

v. -/ 



Hi 



8tt 



(y,+r r ) i(y-y,) f<g 3 

-i(7,~Tl) (7<+7r) N fi, 



2y 



■J 



(K 



(l+0|)[|fr 2 + i(co + s 2 )] 

= (1+7) 

[l + (l+ T )|][(l + 7>3+ 



K =l+7 



0=1+7 



(to— S 2 ) 
(w+Si) 

(co+s 2 ) 

(co — Si) 

m 2 

Mi 



I co J 



X- 



(l-«i)7 (l-fe)? . 



(l— «.)? 



■(.^y •>*!)' -[i + o-w>g 



(44) 



(46) 
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Thus we have been able to express the electrical 
conductivity tensor containing ion-motion, electron- 
electron and electron-proton collisions into a form 
which is similar to that of the Lorentz tensor [Sen 
and Wyller, 1960a]. 

When 7 goes to zero, (45) reduces to a form of 
the electrical conductivity tensor [Wyller, 1961] valid 
for high frequencies (conventional magneto-ionic 
theory). 

The dielectric tensor is derived from the relation 






6 = 1 + ?^ a 

10) 



(47) 



for which we find 

«i=(l— On)— *&u, «n=i(/ii— du)+i(en— en), 

«m=[<hi-i(Cu+«u)]+*l&u-i(/n+4i)], (48) 
where 

«4(l+7) 



«n = 



bn= 



Cir- 



[l + (l+7) |] [(l+7) 2 ^+co 2 ] 

^(l+7) 2 "8 , 

co[l + (l+ 7 )|][(l+ 7 ) 2 , 2 + co 2 ] 

Q)l(o:—S2)a 

(l+ a |)[a^+(co-s 2 ) 2 ] 



dn= 



2 2 
O) a Vi 



en- 



a>(l+a|)[a 2 ^ + (co-S 2 ) 2 ] 

cooCc o + frQfl ? 

[/3 2 ^+(co+s 2 ) 2 ] 



•W) 



"(""*«) 



;l+£sin 2 <ft±(.g 2 sin 4 4>-g 2 cos 2 <ft)^ 
.D+E'sin 2 ^ 



(49) 



A=2e 1 (€ I +em), i?=e m (e I +e UI )+«n, C=2« 1 e II , 
Z>=26i, E=2e ln , 



R=- 



g v _ B sin 2 <j>T (B 2 sin 4 4»-(? 2 cos 2 4>V A 



8. 



Ccos$ 



(50) 



Here n is the real refractive index, ck/co is the ab- 
sorptivity, and <j> is the angle between the direction 
of the magnetic field and the direction of propagation. 

5. Zeroes and Infinities of the Refractive 
Index 

In the formula (49) the complex refractive index 

CK 

n=n—i — satisfies a biquadratic equation 



_ 4 . _, (A+B sin 2 <t>) . A 2 +2AB sin 2 4>+C 2 cos 2 <j> 



(D+E sin 2 4>) 



(D+ E sin 2 <t>) 2 



0. 



(51) 



Hence the condition that n may have a zero root 
is that 

A 2 +2AB sin 2 <t>+C 2 cos 2 ^=0 

provided that D+E sin 2 <£^0. After some al- 
gebraic simplifications the above two expressions 
reduce respectively to 

ei[(ei + €in) 2 + en](ei + €ni sin 2 0)=O (52) 

cj+em sin 2 0^0. (53) 

By virtue of the restriction (53), all the zero roots of 
n will be given by the following three equations, 

(a) €!=0, (b) e I +e I1I —ie I1 =Q J (c) e L +e m +ie n =0. 



fn=- 



u 2 ot3 2 v 2 



(i+isf)^ 



^+(co+S 2 ) 2 ] 



These equations are independent of the propagation 

angle "4>". Assigning the appropriate values to 

e 1 from (48) we find that (a) will hold when a n = l 

and b n = simultaneously. The latter condition 

means that either 7=0, i.e., there is no collision, 

or else that 5 3 =1, which gives an imaginary 

-d , . . ,i i tut n j- v r^™i value to the collision parameter Z = v/os and is 

R * et T?L ng + S% Sen ™ d , W y}t T forn } allsm , [^PO] therefore inadmissible. The first condition with 
we find that the formulae for the complex refractive - =Q ieMs cutoff for n t 

index c/u and the polarization It are unaltered in J 

form except that the constants A, B, C, etc., are 

now functions of e^ e n , and e m defined in (48). X= 

This leads us to a new formulation of the hydro- 
magneto -ionic theory which permits the study of 

low frequency modes such as the Alfven mode and Similar treatment to (b) and (c) yields cutoff 
the retarded magneto-acoustical mode with electron- values of n without collisions at 
electron and electron-proton collisions included. 



1+7 



(54) 



For convenience we recapitulate [Sen and Wyller, 
1960b] 



X- 



q+P(i- 7 D 

1+7 



(55) 



100 



and 



respectively. 



X- 



(1-F)(1+ 7 F) 

1+7 



In (54) and (55), X=c4^ 2 and Y=s 2 /a>. The 
cutoff values are not attained when collisions are 
present. These cutoff values will be illustrated by 
the curves for n against X given at the end of the 
paper. 

From inspection of (51) it is apparent that n will 
have an infinite root when 

D+E sin 2 </>=0. 

This is equivalent to the condition 

6 i+ e in sm2 0=0. (56) 

For the case with collisions, the condition (56) turns 
out to be fairly complicated, but, if we neglect 
collisions, it yields a simple relation analogous to 
that of Hines [1957]. 



It may be noted that the condition (57) for infinite 
root of n is dependent on 0. When = 7r/2, an 
infinite root of n occurs at 



X-- 



(1-F 2 )(1-7 2 F 2 ) 
= (1+t)(1-tF 2 ) ' 



(58) 



This value is somewhat different from that obtained 
by Budden [1961], even if an allowance for his 
approximation is introduced in our computation. 
For the low frequency end of the whistler spectrum 
we consider the case S2»w»§i. In this limit 
the condition (56) simplifies to 



{l-^(l+ T )}cos 2 </,+ / 



1+ 



Wo 



-&\ sin' 



4>=0. 

(59) 



{ 



~!(1+T)} 



cos 2 0- 



+ ^-? \ sin 2 0=0. (57) only when co<a> J. 



If the coefficients of cos 2 and sin 2 </> in (59) are both 
negative, the condition for an infinity will not hold; 
which means that the wave can propagate in all 
directions. The first coefficient will be negative if 
co<o> ; also the second coefficient will be negative 



-M 2 —- These bounds on co will 



St+C4 



Y = 100 




200 300 400 500 
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be helpful for the detection of ions with whistlers 
(co<<co ) which are propagated to a large extent in 
the upper atmospheric regions where our present 
model of a fully ionized plasma will be a close 
approximation [Hines, 1957]. 

6. Numerical Applications 

The real refractive index, n, and the absorption 

CK 

factor. — 9 were evaluated from the relations 

CO 



where from 



v 



r=(n—i—j =M—iN 



n=^M+^M 2 +N* 
V 2 



CK 



■yJZ 



We note that our signs in (49) for the complex 
refractive index are opposite to those of the magneto- 
ionic theory; thus our upper (U) and lower (X) 
signs, correspond respectively to the extraordinary 
and ordinary waves of the magneto-ionic theory. 
We have computed the propagation constants n and 
ck/co as a function of Xior the following representative 



values of: 

(i) Y= 100 (whistler mode) ; 

(ii) Y=-yJl/y (lower hybrid frequency) ; 

(hi) Y=l/y (ion gyroresonance) ; 

and (iv) F=10 4 (hydromagnetic mode). 

We will briefly discuss below the above cases. 

Case (i) Y=100 (whistler mode) 

Longitudinal propagation (<£=#). Figures 1 
and 2. The introduction of ion motion has little 
effect on the birefringence, which is only slightly 
decreased. The important parameter for collisional 
effects seems to be the ratio Z/Y. Even for small 
values of this ratio (0.1 in fig. 2), the cutoff in the 
ordinary ray is removed. The absorption of the 
extraordinary ray is increased with increasing Z=v/oo. 
The reverse is the case with the ordinary ray. 

Transverse propagation (<£=7r/2). Figure 3. L 3 
is not drawn in the figure. It follows very closely 
the vertical axis starting from the value n=l at 
X=0 y with a cutoff at x= 0.995. The effect of ion 
motion is negligible for the ordinary ray. It is quite 
marked for the extraordinary. The birefringence 
of the medium is appreciably decreased. The 
collisions remove the cutoff for both the ordinary 
and extraordinary rays. The inclusion of ion motion 
has no effect on the absorption of the ordinary ray. 
For the extraordinary ray, it at first increases and 
then decreases the absorption. The ion motion 
seems to have appreciable effect for transverse 
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propagation. The effect may well be confined 
within a narrow cone. 

Case (ii) Y=^Tjy (lower hybrid frequency). Figures 
4-7. The behavior of the propagation factors for 
the longitudinal case is as in the magneto-ionic 
theory. But for transverse propagation, collisions 
not only remove the cutoff but reverse the trend of 
the real refractive index curve (figs. 6 and 7) for the 
lower sign (ordinary ray of the magneto-ionic theory) . 
How the introduction of ion motion makes possible 
all directions of propagation for this case has been 
discussed at the end of section 5. 



Case (Hi) Y=l/y (ion gyroresonance) . Figures 
8-10. 

Longitudinal propagation (</>= 0) . We shall show 
next how for longitudinal propagation the ion-gyro- 
resonance introduces a singularity in the lower mode. 
Figure 8 shows that even a small amount of collision 
removes this singularity. The wave can propagate 
but is nevertheless subject to high absorption. With 
increasing collision frequency, the absorption of the 
lower mode decreases, whereas that of the upper 
mode increases (fig. 9). 



104 




10 20 30 40 50 60 70 



Figure 7. 




Figure 8. 



Transverse propagation (4>=ir/2). The singu- 
larity is removed in this case and the wave can 
propagate. The cutoff frequencies for the two 
modes for no collisions is discussed below. The 
propagation factors for a high collision case (Z= 
1000) are shown in figure 10. 

Case (iv) Y=10* (hydromagnetic mode). 

Longitudinal propagation (</>=0). Figures 11 
and 12. These are the retarded magneto -acoustic 
mode (lower sign) and the oblique Alfven mode 
(upper sign) of Denisse and Delcroiz [1961a]. Note 



in figure 12 that the oblique Alfven wave does not 
suffer absorption, but that the retarded magneto- 
acoustic mode does. This property may have 
important application in the solar corona [Oster- 
brook 1961]. 

Transverse propagation (</>=7r/2). Figures 13 
and 14. The lower mode suffers a cutoff at X~l, 
as in the magneto-ionic theory, which is removed 
by collisions. Note that the nonabsorptivity of 
the Alfven mode (upper sign) persists even in 
transverse propagation. 
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We shall now give a brief analysis of the ion sign, we get after a little algebra (for Z=0, i.e.. no 
gyroresonance referred to in case (iii) above. From collisions): 
(49) we have for longitudinal propagation (<£=0): 



W 2 =1 _ (1+7 ) 



(l-Y 1 )X-(l+y)X 2 



(i-r 2 )(i-FD-(i+T)(i-r 1 )X 



n 2 =(A J riC)/D. Taking the minus sign and Z=0 

(no collisions), we have n 2 =n 2 =l — X[(l + Y)~ 1 + 

y(l — Fi) _1 ], where Fi=Si/w. There is therefore a 

singularity in n at the ion gyroresonance (Fi = l). 

It can be shown that a finite Z will remove the When Fi=l (ion-gyroresonance) , n 2 -^2-\-yX. Thus 

singularity. For transverse propagation (</>=7r/2), n 2 -^>2, as X—>0. It can be shown that with collisions 



we have n 2 =(A+B±B)/(D+E). Taking the + 



>1, as X-+0. (See fig. 10.) 
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7. Conclusion 

In the present paper we have given a microscopic 
hydromagneto-ionic theory, to Chapman and 
Cowling's second approximation, of a completely 
ionized hydrogen plasma. As a special case for no 
collisions, i.e., v=o, our (49) reduces to the expres- 
sions derived by Astrom [1951] for the longitudinal 
and transverse components. Furthermore the most 



general dispersion formula derived rigorously from 
Denisse and Delcroix's general conditions [1961b] 
gives the same transverse and longitudinal modes 
of propagation as those obtained from our formula. 
As already stated in the introduction we have further 
agreements with the results obtained by Kantor 
[1963] in the limit of no collisions. All these provide 
an indirect check on the internal consistency of our 
theory. 
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Figure 14. 



We note that the treatment can easily be extended 
to a neutral plasma of any degree of ionization. 
In its present formulation, it should have important 
applications to hydromagnetic wave dissipation in 
the solar corona. The treatment of collisions in 
our hydromagneto-ionic theory is broad enough to 
include any velocity dependence of the collision 
cross sections. For applications to hydromagnetic 
wave dissipation in the terrestrial ionosphere, our 
treatment should be extended to a partially ionized 
gas. Burgers' formalism should be adequate for this 
extended treatment, as is evidenced by Pipkin's 



work [1961] on the d-c conductivity of a partially 
ionized gas. 

The numerical applications that have been given 
here are merely illustrative and by no means exhaus- 
tive of the results that can be obtained from the 
theory. 



It is a pleasure to acknowledge our indebtedness 
to Drs. Pfister, Poeverlein, Kantor, Prasad, and 
Lt. Finn for fruitful discussions, and to Mr. Arnold 
Shickman for the computations. 
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